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It is shown that surface of a liquid consisting of several interpenetrating superfluids moving with different 
velocities becomes unstable at some threshold. We demonstrate that the criterion for the onset of the insta- 
bility changes in the presence of disspative interaction between the surface and the environment. Possible 
physical applications of the surface instability are discussed. 
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1. Introduction. Analogs of the classical Kelvin- 
Helmholtz instability are relevant for various physical 
systems. For instance, recently it was proposed [1] that 
the analog of KH instability for a mixture of two su- 
perfluids can provide a possible mechanism for pulsar 
glitches [2], The interest to KH instabilities has also 
been revived in connection with the observation of the 
instability arising on the interface separating two super- 
fluids, 3 He- A and 3 He-B [3]. The results of this exper- 
iment are in excellent agreement with the criterion for 
the onset of the instability derived in [4] . Free surface 
instability of this class was discussed in [5]. 

However, in most interesting physical situations, 
such as neutron stars and spin Bose-Einstein con- 
densates [6], several interpenetrating superfluids are 
present. Such mixtures of superfluids are known to ex- 
hibit Andreev-Bashkin effect [7] — that is, the mass 
current of each condensate is carried by the velocities of 
the others: 
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where v l and ji are the superfluid velocity and the mass 
current of the i-th component of the condensate. The 
denstiy of kinetic energy of such condensates can be 
written in the following form: 
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In the present Letter we show that the free surface 
of a multi-component condensate exhibiting Andreev- 
Bashkin effect becomes unstable when velocities of its 
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constituents reach certain values. There exist two dif- 
ferent thresholds characterizing the instability, the first 
resulting from the dissipative interaction of the surface 
with the environment, which is typically small, and the 
second corresponding to the classical dynamic instabil- 
ity. When the system is between these thresholds, the 
time of instability development is determined by the 
interaction with the enviroment; after passing through 
the second threshold the characteristic time of instabil- 
ity evolution quickly decreases. In some cases it is pos- 
sible that there would be no intermediate region, that 
is, the instability develops very quickly from its very 
beginning. 

We propose that this surface instability should arise 
earlier than the bulk instability described in [1] . Indeed, 
the critical value of velocity obtained in [1] is compara- 
ble with the sound velocity in the bulk, while the critical 
velocity for arising surface instability is determined by 
the field stabilizing the surface and therefore should be 
much smaller. 

There is another possible application of the surface 
instability arising in mixtures of several interpenetrat- 
ing superfluids: this mechanism could be relevant for 
the instability of a Bose-Einstein condensate with a 
multi-component order parmeter (see [6] and references 
therein). 

The paper is organized as follows. First we derive 
the criterion for the thermodynamic instability, then the 
condition for the emergence of dynamic instability is ob- 
tained. We analyze how the presence of the interaction 
of the surface with the environment affects the insta- 
bility criterion. After that we demonstrate that in the 
'shallow water' limit the physics of the instability re- 
mains the same, but can be described in terms of the 
effective ripplon metric. 
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2. Thermodynamic instability. In this section 
we consider a thermodynamic instability arising on a 
free surface of several interpenetrating superfluids mov- 
ing with different superfluid velocities v\ We suppose 
that the surface in equilibrium coincides with the plane 
z = and superfluid velocities v* are equal to u\ As we 
will show in Section 3, this instability results from the 
interaction between the surface of the liquid and the en- 
vironment. The instability emerges when the energy of 
static perturbations becomes negative in the container 
frame (or, equivalently, in the rest frame of the normal 
components, where normal velocity v„ = 0). 

Let us find the energy of the perturbation with 
wave vector k = (k x ,k y ). First of all, the deforma- 
tion ( = a sin kr of the surface causes the perturbation 
(5v* = V$ 4 of the superfluid velocity fields v\ There- 
fore, the free energy functional for the perturbation is 
given by: 

T= l -Jdxdy [Fe+a((d x 2 + {d y 2 )] + 



frequency of some surface wave acquires positive imag- 
inary part. It appears that the threshold for arising 
of the instability depends substantially on the presence 
of dissipation. Even very small interaction of the su- 
perfluid with the environment (for example, it can be 
friction between the surface of superfluid and the con- 
tainer walls or between the surface and normal com- 
ponent) moves the threshold to another value. Let us 
first suppose that the interaction of the surface with the 
enviroment is absent. 

The Euler equation in our case reads 

E if + E PaWovj, = -V f3P Fp. (8) 

The solution of (8) corresponding to a small-amplitude 
surface wave with frequency ui and wave vector k = 
(k x , k y ) can be chosen in the following form: 



v l x = ul+ik^^-^e^A 1 , 



+ ik y e l ^ kr - ult) e kz A l 



(9) 
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where F is the external field stabilizing the surface, e.g., 
gravitational field, and a is surface tension. Since the 
perturbation is stationary, the continuity equation for 
each component reduces to A$* = 0. Hence $ l has the 
following form: 



= ^expfczcoskr. 



(4) 



We can express constants A 1 via a using the boundary 
conditions 



dr dz 

This gives A 1 = a(u*k), where k = ^. Substituting 
these values of A 1 into the free energy (3), one obtains 
the energy of the perturbation: 

E{k) = ^[F + k 2 a-kJ2pij(^)(^)}- (6) 

i,j 

Thus, the criterion for the onset of the thermodynamic 
instability is given by 



i max pij (u*k) (u J k) = VoF. 



(7) 



3. Dynamic instability. Now we proceed to 
the investigation of the analog of the classical Kelvin- 
Hclmholtz dynamical instability for several interpene- 
trating superfluids. This instability emerges when the 



v l z = kS^-^e kz A\ 



(11) 



p = -Fz + ie i( - kr -^e kz x 

x^p^'-^^Kk)^], (12) 

where A 1 , i = 1, . . . , n are some constants. Apart from 
equation (8), our system must satisfy boundary condi- 
tions: 



p - a ( — + —) C 

\dx 2 dy 2 J 

c i id( d( 
SVz - U d-r = d-f 



(13) 



(14) 



Using equation (13) we can express the surface defor- 
mation ( in terms of constants A 1 : 

( = ie^ r - ut) e kz /(F + k 2 a)x 

-E/'-.-- 1 ' E/v :u ' kLl ' ( 15 ) 

i,3 i,3 

Substituting (15) into (14) and making use of (11), we 
obtain a set of equations for constants A 1 : 

k(F + k 2 a)A l = (w - u 4 k)[^ PjiA j (to - u'k)], 

j,l 

i = 1, ... ,n. (16) 
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It is convinient to rewrite equation (16) in the form 

5^«A''=0, i = l,...,n, (17) 



where 



fa = fc(F + k 2 a)Sij - (w - u*k) ^ p j7 (w - u'k). (18) 

For the system (17) to have a non-trivial solution it is 
neccesary that 



det <j>ij = 0. 



(19) 



It can be shown (see Appendix A) that the condition 
(19) is equivalent to the following equation: 

k(F + k 2 a) = ^p ij (w-u i k)(u;-u J 'k). (20) 



The equation (20) determines the ripplon spectrum. 
The frequency of the ripplon with wave vector directed 
along k acquires imaginary part when 

p]>>>*k)(u^k)-(^>*k) =2 P V¥^, (21) 

i,j \ i,3 / 

where p — J2i j Pij- The equation (21) provides the cri- 
terion for the onset of the instability which is analogous 
to the classical Kelvin-Hclmholtz instability. In terms 
of 'mean velocity': 



jtot _ Pij U * 



(22) 



P P 
equation (21) becomes more transparent: 

]T pa Kk) (u'k) - p(vk) 2 = 2 Vfw. (23) 

As was pointed out in [4], the criterion (23) will 
change if one takes into account dissipative force arising 
when the surface is moving with respect to the container 
walls. This force can be written in the following form: 



F b = -Td t C 



(24) 



In the presence of friction, the equation (20) for the 
spectrum of ripplons modifies as follows 

k(F + k 2 a) - iTuk = Pij(" - u l k)(w - u J k). (25) 



When the supcrfluid velocities reach critical values, 
imaginary part of frequency crosses zero. Then from 
equation (25) it follows that the real part of frequency 



must also cross zero. Consequently, the instability con- 
dition is 



- max ^ p^ (u*k) (u'k) = V oF . 



(26) 



This equation coincides with the condition (7) for the 
thermodynamic instability. From equations (23), (26) 
it can be easily seen that thermodynamic instability 
always arises earlier than the classical dynamic insta- 
bility. In the region between these two thresholds the 
instability develops rather slowly (because the time of 
instability development is proportional to the friction 
parameter T, which is usually relatively small), and af- 
ter passing through the dynamic threshold it becomes 
quickly developing. 

Criteria (23), (26) can be represented in another 
form, which sometimes proves more suitable. In order 
to obtain this form, we note that instability occurs first 
for the ripplons with wave vector fc = \jFja. We fur- 
ther rewrite equation (25) with k = fcok in terms of the 
effective ripplon metric: 



9^K = iTf Q , 
k» = {-iv/k ,k a ), a = x,y. 



(27) 



The effective ripplon metric in equation (27) is given by: 

.,00 
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(28) 



9 a * = 2^8"? - Y.ijfHi^. 

It can be easily seen that criterion (23) for dynamic in- 
stability is equivalent to the condition 



det^,, = oo, 



(29) 



while criterion (26) for thermodynamic instability is 
equivalent to 

5oo = 0. (30) 

4. Case of two interpenetrating superfluids. 

Now let us consider in more detail the case when only 
two interpenetrating superfluids are present. The equa- 
tion (25) for the spectrum of ripplons reads: 



F + k 2 a uj co , « 2 

+ P22 (^ u 2 k) 2 + 2 P12 (^ - u : k)(^ - u 2 k). (31) 

It is convenient to rewrite equation (31) in terms of 
mean velocity (22) and relative velocity u = u 1 — u 2 : 



F + k 2 <r 
k 



- lV k= P 



Q - vk) 2 + a(uk) 2 



(32) 
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where a = (P1P2 — p\i)l P 1 - Thus, the frequency of the 
surface wave with wave vector k satisfies the following 
equation: 



w c , F + k 2 a Flu ~ 

i =vk± V~ — l -pk~ a{uk) ■ 



(33) 



Using (33), one immediately obtains that the first in- 
stability threshold, corresponding to the onset of the 
thermodynamic instability, is determined by the equa- 
tion: 



2-JFo 



max 

k 



a(uk) 2 + (vie) 2 



and the second instability threshold is given by 



2VFa 



au 



(34) 



(35) 



Note that in the absence of Andreev-Bashkin effect re- 
sults (34), (35) coincide with those for the instability of 
the interface separating two one-component superfluids 
[4]. The only difference between these situations con- 
sist in the fact that the stabilizing field F in the case of 
surface instability is proportional to the total density of 
the mixture, while in the case of interface instability it 
is proportional to the difference of the densities of the 
separated liquids. 

As we have already mentioned, the thermodynamic 
instability always occurs before the dynamic one. How- 
ever, under certain circumstances these instabilities can 
arise simultaneously. Indeed, if u is perpendicular to v 
and v < y/au, then equations (34) and (35) coincide. 
Therefore, the instability starts at u 2 = 2v J^ and it 
is 'strong' in the sense that the time of its development 
quickly decreases as u grows. This scenario is analogous 
to the violation of the 'cosmic censorship' principle in 
quantum liquids [8], when the effective ripplon metric 
acquires naked singularity. 

5. Shallow water limit. In the previous sections 
we supposed that the depth of the liquid h is large, so 
that 



/i> i. 



In this section we demonstrate that in the limit of 'shal- 



low water' (i.e. \J F/ah <C 1) the mechanism of the in- 
stability remains the same. Indeed, the equation for the 
ripplon spectrum in this limit has the following form: 

^ / !jy(u;-u i k)(a;-u J 'k) = Fhk 2 +<rhk 4 -ir{k)uj. (37) 



It can be rewritten in terms of the effective ripplon met- 
rics: 

g^k^K = iT(k)uj - <rhk A , (38) 



fc p = (-^,fc a ), a = x,y. 



g m = -p, g 0a = -j totc 



g a e = Fh8 af3 



(39) 



(40) 



(41) 



Note that in contrast to the metric (28) the metric (40), 
(41) is valid for all long-wavelength ripplons. 

Again, there would be two thresholds for the onset of 
the instability, which are completely analogous to (23), 
(26), the only difference being that the ripplons which 
are responsible for the instability have wave vectors 
close to zero. The thermodynamic instability thresh- 
old corresponds to the condition (30) goo = 0, which 
ensures the appearance of the horizon, as was recently 
proposed by Schutzhold and Unruh [9]; the condition 
for the onset of the dynamic instability is detg^,, = oo, 
in agreement with [8]. 

6. Conclusion. To conclude, we have found that 
the free surface of several interpenetrating superflu- 
ids, moving with different velocities, becomes unsta- 
ble at some threshold. The criterion for such instabil- 
ity contains the off-diagonal densities pij , characterising 
Andreev-Bashkin effect. This kind of instability could 
provide a triggering mechanism for pulsar glitches. It 
also probably takes place in atomic Bose-Einstein con- 
densates. 
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Appendix A. Here we present a calculation of the 
determinant of the matrix (18). It proves useful to in- 
troduce the following notations: 



a" 



u*k, 



(5 l =5>^-u 4 k), 

i 

7 = fc(F + fc 2 cr). 
Then matrix (18) takes the form 

/ 7 -a 1 /3 1 -a 1 /? 2 ••• ~ 
-a 2 /? 1 1 -a 2 l3 2 ... 



(42) 

(43) 
(44) 



y -a"/? 1 -a n (3 2 



-a 2 (3 n 
~f-a n (3 n j 



(45) 
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The determinant of matrix fcj can be written as follows: 



det <f> = TT ct x 

i=l 

/ 7/a 1 -/? 1 -/3 2 ... -/3" \ 
-P 1 7/a 2 -/? 2 ... -/?" 



x det 



v 



-/3 2 ... j/a n -/3 n / 



(46) 



In order to calculalate the determinant of the matrix on 
the r.h.s. of (46), let us subtract its (n — l)th line from 
its nth line, its (n — 2)th line from its (n— l)th line, etc. 
This gives 



D n 



det 



— 7/a 1 7/a 2 











V 



(47) 



-7/a"- 1 7/0;" / 

One can easily obtain a reccurent formula for D n : 



i=i 



Thus, 



(48) 



(49) 



det(/. = 7™- 7 ™- 1 ^a l /3 J , 

i=l 

and it equals zero when 

n 

7 = 5> 2 /? 1 . (50) 



i=l 
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